CERN-TH/99-134 
hep-ph/9905380 



Low-energy neutrino-photon inelastic interactions 



JOAQUIM MATIAS * 

TH-Division, CERN, Geneva 23, Switzerland 

a\ 

G\ 

>v Abstract 

5 

qq ■ The computation of the polarized amplitudes and cross section of the pro- 

cesses 71/ — ► 77^, 77 — > 71^ and i/p — > 777 is described. We used an 
effective lagrangian approach for energies below the threshold for e + e~ pair 

<^> , production and the complete computation at higher energies for application 

OO ■ 

. in supernova dynamics. Leading contributions of physics beyond the SM are 



in 
o 

On 
On 

a: 

Oh! 



X 



also commented. 



May 1999 



*Talk given at XXXIV* h Rencontres de Moriond: Electroweak Interactions and Unified Theories, 
Les Arcs, France, 13-20 Mar 1999. 



1 



I. INTRODUCTION 



Low-energy neutrino-photon interactions are of potential interest in astrophysics, since 
they could affect the mechanism of stars to lose energy and hence the study of the stellar 
evolution, in particular, supernova dynamics. 

However the cross sections of the 2 — > 2 processes (77 — > vv, jv — > ^v and vv — > 77) are 
too small. The reason for their strong suppression is the prohibition of the coupling of two 
photons to a J = 1 state of any parity. This is simply Yang's theorem [Ij] . As a consequence 
the amplitudes of the 2 — > 2 processes are zero to order Gf- For instance [fj], 

A xx {iv -> 71/) = 1 + - log — £ cos —f , (1) 



2tt V 3 to / 

where A, A' are the photon helicities and f xx ' is a certain function of s and t. The important 
point to notice here is that the scale of the process is given by M w , which strongly suppresses 
it at low energies. However there is a way to bypass this suppression: if one couples three 
photons instead of two the theorem does not apply anymore. This implies on the one hand, 
that we will pay an extra a in the cross section. But, on the other, that it is possible that 
the scale of the process is no longer Mw but some other light scale that could give a large 
enhancement. 

The processes will then be 

77 — > vvj, 'jv — > 77Z/, and vv — > 777 (2) 

The first may affect the stellar energy loss mechanism and the other two may reduce the 
mean free path of a neutrino inside a supernova core and in principle they could act as a cut 
off of high energy photons or neutrinos by the background. The only existing computation in 
the literature of one of these complicated processes disagrees || with a recent computation 
by Dicus and Repko |4[|. Our aim in || was to settle down this disagreement and give 
an explicit derivation of this effective lagrangian that could be useful for applications in a 
different context. 

Numerically, it was found that the cross section of the inelastic 5-leg processes are be- 
tween 9 and 13 orders of magnitude larger than the 2^2 corresponding processes for 
center of mass energies u between 0.2 and 2 times the mass of the electron m e . 

However if one needs to go to higher energies, for instance to energies above the threshold 
for e + e _ pair production, an exact calculation of the processes in Eq. (0) was important so 
as to definitively assess their role in astrophysics and the range of validity of the effective 
theory. 

In this talk we will comment on both approaches: the effective lagrangian and the direct 
computation in the SM. We will also briefly comment on how physics beyond the SM could 
affect them. 
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II. EFFECTIVE LAGRANGIAN APPROACH 

Three main observations will allow us to obtain the leading contribution at low energies 
(of the order of m e ) to the cross section of the processes in Eq. (|2|): 

(i) First, the lack of the 1/M^ suppression in the 5-leg processes (Yang's theorem does 
not apply there) suggest that we concentrate on those diagrams with a lighter particle 
inside the loop, for instance m e . 

(ii) Second, the processes in Eq. (|2|) are finite, no counterterm is needed. This observation 
helps in the search for the leading diagrams by defining a hierarchy of diagrams at low 
energies. Giving a certain topology it becomes easier to find its large m e limit. These 
two observations imply that the leading diagrams contributing to these processes are 
those given in Fig. 1. Any other diagram or topology will include other particles 
inside the loop, different from the electron, and they will be automatically suppressed 
by inverse powers of its mass. 




Figure 1. SM leading diagrams contributing to five-leg photon-neutrino processes: a) Type 

A diagrams b) Type B diagrams 

(Hi) So far we have found the leading diagrams. A third observation will allow us to obtain 
its low-energy expansion in an elegant way: the subset of diagrams depicted in Fig. 1 
resembles a photon-photon scattering process (see Fig. 2), one of the photons being 
off-shell. 



P2, e 2 




Pi, e 4 

Figure. 2: Four-photon interaction: Type C diagrams 
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Before expressing this "resemblance" in mathematical terms we should clarify why it is 
interesting to make this link. The reason is as follows: we know that the SM is able to de- 
scribe both the photon-photon scattering process and our inelastic photon-neutrino process 
which we want to compute. Moreover, we know that the low-energy limit of the photon- 
photon scattering process at one loop is given by the Euler-Heisenberg lagrangian || . Hence 
if we are able to find a connection between the diagrams of the two processes, we will auto- 
matically be able to establish a link of the effective lagrangian of our process in terms of the 
Euler-Heisenberg lagrangian of the photon-photon scattering without having to compute a 
single one-loop diagram. 

In order to establish the link between those two processes, we should be able to answer 
two questions: 

(a) Can type B diagrams be reduced to type A diagrams? 

(b) Is it possible to relate type A diagrams, which contain a vertex Zee with an axial part 
Wlfj,{ v e + <V75)/(2c0), with type C diagrams whose corresponding vertex igsej^ does 
not? 

It was shown in ]5] that it is possible to give a positive answer to both questions. The 
keypoints of the proof are two: on the one hand, the gauge boson propagators should be 
expanded in the large Mw, Mz limit. This is justified since we are working at very low 
energies. On the other hand, we should find out the combination of the amplitudes of type 
A diagrams with different polarization indices in such a way that the axial part cancels. 
Indeed, Gell-Mann gave an indirect proof in [/]], using charge-conjugation arguments, that 
such combinations should exist. 

In conclusion, at leading order in it was found in Q that the following set of four 

diagrams (from a total of 12) of type A and type B diagrams is proportional to the same 
integral, called Li (see || for definitions): 

5 3 i 
^123 + ^321 + n l23 + -D321 — 2 — ^ e > >M ~KzC^~ ' 

Similar results are obtained for the other two groups of four diagrams, the only difference 
being a trivial change of momenta and indices inside . 

At this point, the correspondence with the four-photon scattering is evident. In fact, by 
fixing the fourth photon leg and calling C^ 7 the corresponding amplitude in Fig. 2, one 
finds a contribution proportional to the same integral L\\ 

cfls + cgf = 2g 4 sW&, A 4 )L^ 7 , (4) 

and similar results for the two remaining combinations Cf 3 ^ 7 + C^j 7 and C^ 7 + C^ 7 . 



4 



Up to now, we have proved that both processes are governed by the same integral, so 
both should have the same momentum dependence. Moreover, we have said that at low 
energies photon-photon scattering is governed by the Euler-Heisenberg lagrangian || . This 
automatically tells us which are the only operators that can be generated at low energies. 
We have no freedom either in the structure nor in the relative coefficients between the 
operators. Equations (|3|) and M) suggest to define a new gauge field up to the neutrino 
current A v = — 75)^ = 2T U , with field strength F^ v . Indeed this definition is not 

unique, so we should add a constant in front of the effective lagrangian that still has to be 
fixed. The lagrangian then reads 

C 
180 

In order to fix this constant, we can use the following ratio of amplitudes 

ylSM /Jeff 

lim :2±L = (6) 

larger Af ' K) 

where P is our process and 47 is the photon scattering. Then we obtain || 

g 5 s 3 w (l + v e ) = 2G F a^ 2 {l + v e ) 

32n 2 mjM 2 v V^mj ' 1 ' 



eff 



5 (F, U F^) (f Xp F x p) - UF^F» x F Xp F^} . (5) 



At this point some remarks are in order. Notice that the suppression factor for the 5-leg 
processes is 1/Myy, to be compared with the suppression of the processes with one 

photon less. Moreover the scale is not given by My/ but m e , giving rise to an important 
enhancement. Finally, once we have established the dictionary that translates between the 
two processes, C — > ^ and T M — > 2e A1 (P4, A4), one can check each polarization amplitude 
with the corresponding one of the photon-photon scattering after putting all the photon 
legs on-shell. This provides a non-trivial check of our computation. 

In we obtained, using this effective lagrangian, the eight polarized differential cross 
sections for the processes — > 7Z/7 and 77 — > uuj. Once added we compared our result with 
the unpolarized differential cross section and the total cross section given in |Q. We found 
perfect agreement and hence disagreement with the ones obtained by Hieu and Shabalin [El . 



III. DIRECT COMPUTATION IN THE SM AND BEYOND 

For energies above the e + e~ threshold it is required to perform a direct computation in 
the SM HI! . In order to compute the processes in Eq. (§) directly in the SM we need to face 
a multileg computation. The traditional methods, such as tensorial decomposition [Tl| , suffer 
from large numerical instabilities due to the large proliferation of terms and, in particular, 
due to the appearance of gram determinants A = detkikj that vanish in collinear regions 
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of the phase space where the cross section is, indeed, well defined. We have used a new 
method, specially suited for this type of analysis [11]. This is a modified version of the one 
of Campbell, based on grouping the coefficients in sets with a well defined A — > limit. 
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Figure 3: — > 77^ cross section in fb as a function of u>/m e . 



The keypoints of the method are mainly two: 

(1) The use of a specific representation for the polarization vectors. 

(2) By means of 7-algebra and spinor manipulations, we can reconstruct in the numera- 
tors the structures that appear in the denominators rather than making a tensorial 
decomposition. In that way all the results are expressed in terms of scalar functions 
with 3 and 4 denominators, rank-1 integrals with 3 and 4 denominators, rank-2 inte- 
grals with 3 denominators, and rank-3 functions with 3 denominators. This already 
provides an important simplification with respect to the standard decomposition, in 
that the computation of tensors such as 



r-ppv, pvp\ pupa 



d n q 



q p q u. 



D D^D 2 D {23) 



is completely avoided. 



(8) 



Such 



We made a large use of the Kahane-Chisholm manipulations over 7 matrices JT2[ 
identities are strictly four-dimensional, while we are, at the same time, using dimensional 
regularization. Our solution || is splitting, before any trace manipulation, the n-dimensional 
integration momentum appearing in the traces as [11] q —>■ q + q, where q and q are the 
four- dimensional and e-dimensional components (e = n — 4), respectively, so that q ■ q — 0. 
The net effects are, on the one hand, that the 7 algebra can then be safely performed in 
four dimensions and, on the other hand, that a set of extra integrals containing powers of 
q 2 in the numerator arise, but they are straightforward to compute. 
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As an example, we have plotted in Fig. 3 the result of the direct computation for the 
cross section of the second process || in Eq. (||]). The result is in full agreement with those 
reported in ||. From the plot it is also clear that the effective theory is valid only when 
uo < 2m e , as expected. Since the exact formulae are too involved to be given explicitly, 
we followed two approaches. We first tried to extend the validity of the effective theory by 
computing the next-to-leading term. While this was a nice and completely independent check 
of the results of the effective theory, we found that it is not sufficient to enlarge the range of 
validity. A second approach was to fit our curves. For instance, o[^v — > 77^) = a c?i [^v — > 
771/) x r - 2 - 76046 x exp [2.13317 - 2.12629 log 2 (r) +0.406718 log 3 (r) -0.029852 log 4 (r)] , where 
the effective cross section <r eff is given in Eq. (26) of ||. The range of validity is now between 
1.7 < r = oj/m e < 100. These fits are useful when these results are used for simulations 
in supernovae dynamics. Finally, in [IT3J we studied the leading contribution to the process 

— > 77Z/ in supersymmetry with R-parity breaking and in left-right symmetric models 
(LRSM). The observations (i) and (ii) in section 2, also apply here. The leading additional 
contribution will come in a LRSM through the substitution of the W ± and Z propagators 
by the corresponding W'^ and Z', while in the supersymmetric case it will come through 
the substitution of the W ± by a slepton. Indeed, it is interesting to notice that in the second 
case the lepton number is no longer conserved, i.e. transitions into different neutrino species 
are allowed. In particular, we found that in this process a muon-neutrino prefers to convert 
into a tau-neutrino rather than an electron-neutrino, in agreement with SuperKamiokande 
results. However, our processes hold for energies much below the energy of the atmospheric 
neutrinos (> 1 GeV) and the cross sections are still too small. We found that the cross 
sections in the fip MSSM are enhanced by a factor of the order of few 10% at best, relative 
to the SM cross sections, while the correction in LRSM models is negligible. 



IV. ASTROPHYSICAL IMPLICATIONS 

While the process 77 — > ^vv provides an energy loss mechanism for stellar process and, 
in particular, it could be important in the cooling of neutron stars, the other two processes 
1/7 — > ^77 and vv — > 777 affect the mean free path of a neutrino inside a supernova ans 
should be included in the supernova codes. 



In |L4| it was found, using a Monte Carlo and the results of the effective theory [IH that 
for a large range of values of the temperature (T) and chemical potential {p) the mean free 
path was less than the size of a supernova core (10 6 cm). This result remains valid when 
using the exact computation for certain values of \i and T and at energies not to far from 
the e + e~ pair production. However, being the exacts results available it would be of extreme 
interest to find out the precise values. 

Also in it was predicted using the data of the supernova SN1987A that the exponent 
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of the energy dependence in the cross section a(jv> — > 77^) oc u 1 should drop out from 
7 = 10 to less than 8.4 for uj a few MeV. Using our results we confirmed this prediction and 
indeed, we found that the exponent drops around 3 for energies between 1 to 10 MeV. 

Finally, concerning the possible cosmological implications it was suggested in M that the 
process 71/ — > jjv might be relevant for some cosmological considerations if the decoupling 
temperature |TJ| in the exact case is found to be low enough. With the exacts results we 
found m that the temperature is too high to be of any relevance in cosmology. 
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